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PLASMA HEATING DURING A CORONAL MASS EJECTION OBSERVED BY SOHO 
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ABSTRACT 

We perform a time-dependent ionization analysis to constrain plasma heating requirements during 
a fast partial halo coronal mass ejection (CME) observed on 2000 June 28 by the Ultraviolet Coron- 
agraph Spectrometer (UVCS) aboard the Solar and Heliospheric Observatory (SOHO). We use two 
methods to derive densities from the UVCS measurements, including a density sensitive O V line ratio 
at 1213.85 and 1218.35 A, and radiative pumping of the O VI AA1032,1038 doublet by chromospheric 
emission lines. The most strongly constrained feature shows cumulative plasma heating comparable 
to or greater than the kinetic energy, while features observed earlier during the event show cumulative 
plasma heating of order or less than the kinetic energy. SOHO Michelson Doppler Imager (MDI) 
observations are used to estimate the active region magnetic energy. We consider candidate plasma 
heating mechanisms and provide constraints when possible. Because this CME was associated with a 
relatively weak flare, the contribution by flare energy (e.g., through thermal conduction or energetic 
particles) is probably small; however, the flare may have been partially behind the limb. Wave heating 
by photospheric motions requires heating rates significantly larger than those previously inferred for 
coronal holes, but the eruption itself could drive waves which heat the plasma. Heating by small-scale 
reconnection in the flux rope or by the CME current sheet is not significantly constrained. UVCS line 
widths suggest that turbulence must be replenished continually and dissipated on time scales shorter 
than the propagation time in order to be an intermediate step in CME heating. 

Subject headings: Sun: activity — Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: UV 
radiation — techniques: spectroscopic — magnetic reconnection 



1. INTRODUCTION 

The understanding of astrophysical phenomena gener- 
ally begins with the energy budget. The energy bud- 
gets of coronal mass ejections (CMEs) include contribu- 
tions from magnetic energy, bulk kinetic energy, ioniza- 
tion energy, gravitational potential energy, thermal en - 
ergy, and energetic particles (e.g., lEmslie et al.l [2005). 
Energy can be lost to the system through radiation and 
thermal conduction. White light observations from in- 
struments such as the Large Angle Spect roscopic Coro- 
nagraph (LASCO; iBrueckner et al.l fl995h on board the 
Solar and Heliospheric Observatory (SOHO) allow a 
straightforward determination of the c ontributions by 
kinet i c and gravitational energy (e.g.. IVourlidas et al.l 
I2000L 120101 : iSubramanian fe Vourlidasl I2007D . The mag- 
netic energy is thought to be the largest component of 
the CME energy budget but is difficult to constrain via 
remote sensing because of the lack of good magnetic field 
diagnostics, especially for transient events such as CMEs. 
However, the total and free magnetic energy of precur- 
sor active regi ons can be estimate d using vector magne- 
tograms (e.g.. IMetcalf et aLl 12005). line-of-sight magne- 
tograms from the Michelson Doppler Imager (MDI) on 
SOHO, and empirical rela tionships between X-ray lumi- 
nosity and magnetic flux (Fis her et al.lll998[ ). The mag- 
netic field structure of interplanetary coronal mass ejec- 
tions (IC MEs) can be investigated using in situ measure- 
ments (see lZurbuchen fc Richardsonl2006l and references 
therein) by spacecraft such as the Advanced Composition 
Explorer (ACE). 

Several different lines of evidence now suggest that the 
thermal energy input into CMEs is comparable to the 
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kinetic energy of the ejected plasma. First, observations 
between 1.5 and 3.5 solar radii by the Ultraviolet Cor on- 
agraph Spectrometer rUVCS: lKohl et al.l[l995l 120061 ) on 
board SOHO have been analyzed using a ti me- dependent 
ionization code for fo ur prior ev ents (jAkmal et al. 2001; 
iCiaravella et al.ll2001l : iLee et al.ll200l lLandi et al]l2010f ) . 
The thermal energy input was constrained to be com- 
parable to or greater than the kinetic energy for several 
features during each of these events. This method is used 
in this paper and further described in Section [5] Second, 
in situ measureme nts of ICMEs typi c ally show high ion 
charge states (e.g., ILepri et al.ll2001l : iLvnch et al.l 120031 : 
ILepri fc Zurbuche n 2004), althoug h low charge state 
plasm a is sometimes observed (e.g., [Lepri & Zurb uche"nl 
2010). Detailed models of ICMEs by iRakowski et al.l 
(200?]) showed that continual heating was required out 
to several solar radii to be consistent with the observed 
charge states of both iron and oxygen at 1 AU, with to- 
tal heating requirements again fo und to be comparable t o 
the C ME kinetic energy (see a l solGruesbeck et aD2011f ). 
Third, iFilippov fc Koutchmyl (j2002f ) presented observa- 
tions of a rising prominence by the Transition Region 
and Coronal Explorer (TRACE). During the eruption, 
the feature seen at 171 A suddenly changed from absorp- 
tion to emission by the pro minence, indicati ng significant 
and rapid heating. Fourth, iLiu et all ([2006) investigated 
in situ measurements between 0.3 and 20 AU and found 
that dissipation of turbulence can explain the observed 
heating rates. The turbulence generation mechanism is 
not understood in the inner heliosphere but heating by 
pickup ions contribute to ICME heating in the outer he- 
liosphere. 

Several analytical m odels explore th e energ etics of ex- 
panding flux ropes. iKumar fc Rustl (jl996l . hereafter, 
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KR) assume global conservation of mass, magnetic flux, 
and helicity for a self-similarly expanding force-free flux 
rope. The magnetic energy is found to decrease monoton- 
ically during the process of expansion. Some of the lost 
magnetic energy goes into overcoming solar gravity and 
increasing the bulk kinetic energy, but a large fraction is 
presumed to go in to heating the p lasma within the ex- 
panding flux rope. I Wang et al.l (j2009[ ) develop a different 
model of a self-similarly expanding flux rope which re- 
laxes several assumptions made by KR and implicitly in- 
cludes some effects associated with the solar wind. Their 
inference of a CME polytropic index of | for one event 
indicates continued heating duri ng flux rope expansion. 
iLvutikov fc Gourgouliatosl ([2010D model CME flux ropes 
as expanding force- free spheromaks. By allowing for fi- 
nite dissipation during expansion and assuming a large 
anoma lous resistivity ( e .g., du e to wave-particle interac- 
tions), iRakowski et al.l ([201 If ) show that this model can 
reproduce charge states observed within the flux rope, 
but not the higher charge states such as Fe 16+ in the 
plasma trailing the flux rope. 

However, none of these models specify the physical 
mechanisms responsible for plasma heating. The most 
likely candidate is dissipation of magnetic energy. There 
are several candidate mechanisms, which are discussed 
in detail in Section [51 These include upflow from the 
CME current sheet, kinking of the CME flux rope, small- 
scale reconnection or tearing behavior within the CME 
ejecta, waves driven by either photospheric motions or 
by the eruption itself, thermal conduction, and depo- 
sition of energy by e nergetic particles. Additionally, 
iFilippov fc Koutchmvl ((2002) argue that colliding flows 
along flux tubes can heat prominence plasma to coronal 
temperatures in upward concave regions of flux tubes. 

In this paper we use a time-dependent ionization anal- 
ysis to constrain the thermal energy content and plasma 
heating of a fast CME observed by SOHO on 2000 June 
28. This technique has previously been used to in- 
vestigate heating during four slow CMEs observe d by 
SOHO/XJYCS (see Table [Q). lAkmal et al.l (|200l an- 
alyzed a CME on 1999 April 23 and found that the 
cumulative heating energy was comparable to the ki- 
netic and gravitational potential energy of this event. 
There was a core of cool plasma radiating in C ill which 
was surrounded by ho tter material radiating in O VI. 
iCiaravella et al.l ([200 ll ) analyzed a CME on 1997 Decem- 
ber 2 and found that gradual heat release mechanisms 
are probably more appropriate than mechanisms where 
heatin g is concen trated during the early evolution of a 
CME. ILee et al.l (j2009[ ) investigated a CME observed on 
2001 December 13 and constrained the cumulative heat- 
ing energy for the bright knots observed in O VI t o be 
greater than the kinetic energy. lLandi et"aTI (|2010f ) an- 
alyzed the 'Cartwheel CME' observed by SOHO, Hin- 
ode, and STEREO on 2008 April 9 and again found 
that the cumulative heating energy was constrained to be 
great er than the kinet ic energy. Both ILee et al.l ((2009) 
and lLandi et al.l (|2010f ) show that thermal conduction is 
probably too slow to sufficiently heat the CME plasma 
during the ea rly evolution of the flux rope. Using differ- 
ent methods, Be mporad et al.l ([20071 ) investigated a slow 
flareless CME observed by UVCS and the Large Angle 
and Spectrometric Coronagraph Experiment (LASCO) 
on 2000 January 31 and found that the core and leading 



TABLE 1 

Census of CME ene rgy budget papers us ing the technique of 
IAkmal etTlI (pOOll) 



Event 


Vpos (km/s) 


Mass (g) 


Reference 


1997 Dec 12 


211 


2.9 x 10 14 


Ciaravella et al. (2001) 


1999 Apr 23 


523 


2.7 x 10 15 


Akmal et al. (2001) 


2000 Jun 28 a 


1198 


7.3 x 10 15 


Murphy et al. (2011) 


2001 Dec 13 a 


864 




Lee et al. (2009) 


2008 Apr 9 


650 


1.4 x 10 15 


Landi et al. (2010) 



Partial halo CME; uncertain mass. 



edge of the CME were hotter than the ambient corona. 

The 2000 June 28 CME analyzed in this paper has 
been previously investigated i n multiple separate works. 
iRavmond fc Ciara vella (2004) used radiative pumping of 
the O VI AA1032,1038 doublet observed by UVCS to 
find number densities ranging from ~1.28 x 10 6 -4 x 10 7 
cm -3 . This method is described further in Section [5] and 
is used in this work to find densities for sev eral features 
along the UVCS slit. ICiaravella et~aT1 (|2005h used UVCS 
and LASCO observations to identify the leading edge 
of this CME as a fast mode shock front. Ma ricic et al.l 
(2006) studied the kinematics of the prominence and 
leading edge of this CME with LASCO and the Mark 
IV (MK4) coronagraph at the Mauna Loa Solar Obser- 
vatory (MLSO). This CME has also be en identified as 
a solar energetic particle (SEP) event ([Ho et al.l 120031 : 
iWang et al.1l2006D . 

The organization of this paper is as follows. In Section 
[2j we discuss the methods we use to find densities using 
UVCS spectra. Observations of the 2000 June 28 CME 
are described in Section[3] Section [4] identifies UVCS fea- 
tures that have good density diagnostics. The method for 
our time-dependendent ionization analysis of this event 
is described in Section A discussion of the results of 
this analysis is presented in Section [6l The magnetic en- 
ergy of the precursor active region is estimated in Section 
[71 Constraints on candidate heating mechanisms are dis- 
cussed in Section [U Section [9] contains a summary and 
conclusions. 

2. DENSITY DIAGNOSTICS 

A key component of the time-dependent ionization 
analysis performed in this paper is the electron number 
density at the location in the corona observed by UVCS. 
Knowledge of the UVCS density allows us to exclude a 
significant portion of parameter space when using the 
method described in Section [5] In this section we discuss 
our two primary density diagnostics: a clas sical density 
sensitive O V line ratio ([ Akmal et al.l l2001f ) and radia- 
tive pumping of the O VI doub let by chromospheric lines 
([Raymond fc Ciaravellall2004fl . 

2.1. The density-sensitive AA1213, 1218 O V line ratio 

The ratio of the [O v] forbidden line at 1213.85 A to the 
O v] intercombination line at 1218.35 A is a useful diag- 
nostic for coronal densities of ~10 6 cm -3 ([Akmal et al.l 
|2001[ ) and is based on well-understood atomic physics. 
At high densities the forbidden line is weak because of 
collisional deexcitation. These lines occur at equilib- 
rium temperatures of ~2 x 10 5 K. We use the CHIANTI 
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Fig. 1. The ratio of the [O v] line at 1213.85 A to the O v] line 
at 1218.67 A as a function of electron density at 2.24 x 10 5 K. 



database (jDere et al.l Il997i H009), which includes pro- 
ton collisions. The intensity ratio of these two lines, 
^1213/^1218? is shown as a function of electron number 
density in Figure HJ 

The [O v] and O v] lines straddle Lya at 1215.67 A. 
Bright Lya emission can obscure weak [O v] emission 
to make it difficult to get an accurate line strength for 
O v]. Moreover, bec ause of the construction of UVCS 
(described in Section ET4|) , the secondary channel [O v] 
line sometimes appears at the same location on the detec- 
tor as the primary channel N ill A991.58 A line. There- 
fore, this density diagnostic is most useful at certain 
instrument-dependent line-of-sight velocities and when 
Lya emission is relatively weak, but has the advantage 
that it is weakly dependent on temperature. 

2.2. Radiative pumping of VI 

The second density diagnostic is the intensity ratio of 
O VI A1031.91 to O VI A1037.61. When collisional exci- 
tation dominates, the ratio will be 2:1. Departures from 
this ratio are due to radiative pumping of O VI A1037 by 
C II AA1036.3, 1037.0 near velocities of 172 and 371 km 
s _1 , O VI A1037 by O vi A1032 near velocities of 1650 
km s ~ 1 ,orO VI A1032 by Ly0 near velocities of 1810 km 
s" 1 (jNoci et al.lll987URavmond fc Ciaravellall2004f ). 

The electron number densities are given as follows. For 
pumping of O VI A1037 .61 by C H A1036 .3 near velocities 
of 371 km s _1 (see also lNoci et al.|[l9870 . the ratio is less 
than 2:1 and the number density is given by 



ctio37WC u)W 2 - R 

91037 2i? 



(1) 



For pumping of O VI A1037.61 by O VI A1031.91 near 
velocities of 1650 km s _1 , the ratio is less than 2:1 and 
the number density is given by 



^i037/disk(1032)^2-i? 

#1037 2i? 



(2) 



For pumping of O VI A 103 1.91 by Ly/3 near velocities of 
1810 km s _1 , the ratio is greater than 2:1 and the number 
density is given by 



<ri032/disk(Ly/3)W 2 



(3) 



<2i032 R-2 
In the above relations, <Jio32 and (J1037 are the effective 



scatting cross sections, /disk is the solar disk intensity of 
each of the lines, W = 2tt (l — y/1 — r -2 ) is the dilution 
factor for a distance r from Sun center, #1032 and #1037 
are the collisional excitation rate coefficients, and R = 
^1032/^1038 is the intensity ratio of the two O VI lines. 

This density diagnostic requires more assumptions 
than the O V density sensitive line ratio. The caveats 
of this method are discussed by iRavmond fc Ciaravellal 
(2004) and include that the solar disk intensity of each 
of the illuminating lines may be enhanced early during 
flares ([Raymond et al.l 120071 : Uohnson et al.l 120111 ) , mul- 
tiple components to the plasma can exist at different 
speeds along the line of sight, and the collisional excita- 
tion rate depends relatively strongly on the temperature 
of the plasma. In addition, weak C II AA1036.3, 1037.0 
lines complicate the process of finding the O VI line ra- 
tio from observational data. We mitigate the effects of 
these C II lines by fitting Gaussian line profiles to Ly/3 
and the O VI doublet and enforcing that both lines in 
the O VI doublet have the same width. In Section [4] we 
compare the densities derived through both the O V and 
O VI diagnostics for one particular feature (blob E) and 
find that the results are consistent to within expected 
systematic uncertainties. 

3. DATA SET 

Observations of the CME on 2000 June 28 are available 
from the EIT, LASCO, and UVCS instruments on board 
SOHO, and the MK4 coronagraph at MLSO. The C class 
flare associated with this CME was detected by GOES. 
Pre-CME and post-CME observations of this event were 
made by Yohkoh/SXT and SOHO/MBl. The observa- 
tions for this event ha ve b een previously descr ibed by 
ICiaravella et al.l (|2005D and iMaricic et al.l (|2006f h which 
we summarize below. 

3.1. SORO /EIT observations 

S OHO /EIT images in the 195 A band were taken with 
a ~12 minute cadence during most of this event. Be- 
tween 19:00:14 and 19:19:44 UT, EIT images at 171 A, 
284 A, 195 A, and 304 A were taken in order with a 
six minute cadence. Several of these observations are 
sho wn in Figure |2] fr om 18:00 UT until 19:36 UT (see 
also ICiaravella et al.ll2005l Figure 1). The series of 195 
A observations shows a dark arcade near the northwest 
limb starting to rise around 18:00:10 UT. That it is dark 
at 195 A indicates the absence of plasma at ~1 MK in 
the rising structure. Between 18:36:10 and 18:48:12 UT, 
this feature rises rapidly and is not apparent in the next 
195 A observation at 19:13:48 UT. The 304 A image 
at 19:19:44 UT shows bright strands indicative of He II 
emission, suggesting that the rising filament was rela- 
tively cool. The 195 A images at 19:26:04 and 19:36:24 
UT show a thin strand of hot plasma that outlines the 
He II arch in the 304 A image at 19:19:44 UT, indicat- 
ing the presence of some plasma around 1 MK associated 
with the strand. These strands are probably associated 
with the feature identified in Section [4] as blob F. 

3.2. MLSO/MK4 observations 

MLSO/MK4 measures the polarization brightness 
(pB) of the corona between ~1.1 and ~2.8 solar radii 
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Fig. 2. SOHO/mT observations of the 2000 June 28 CME, 
rescaled to emphasize faint structure. 

using radial scans with a cadence of about three min- 
utes. MK4 observations were performed on the day of 
the event between 16:56 and 20:00 UT. A time sequence 
of the eruption is presented in Figure [3l 

The first clear sign of the CME was at 18:46 UT when a 
rising bright arch entered the field of view (see Figure [3^). 
We identify this arch as the erupting filament observed 
by EIT, most notably at 18:48 UT (see FigureEj). Conse- 
quently, this filament is observed simultaneously by EIT 
and MK4. T he kinematics of this rising prominence were 
presented by Maricic eFall (|2006h . The prominence was 
accelerated most strongly between ~18:40 and ~19:00 
UT at up to ~0.5 km s -2 . The upper parts of this arch 
begin to fade around 19:01 UT (Figure [3b) and are not 




Fig. 3. — MLSO/MK4 polarization brightness observations of 
the 2000 June 28 CME, rescaled to emphasize faint structure. Es- 
timated positions of the features observed by UVCS and discussed 
in Section |4]are labeled. The position of the precursor active region 
is denoted by the red plus sign. 

apparent by 19:16 UT (Figure [St) . 

After about 18:49 UT, the rising loop developed a 
twisted or helical structure. This behavior is consistent 
with the onset and nonlinear growth of a long wave- 
length kink instability, and is apparent in the observa- 
tion at 19:01 UT (Figure [3b). A bright stream of ejecta 
directed towards the southwest is observed to miss the 
UVCS slit entirely (Figure [3b-d). A faint feature rising 
at Vpos ~250 km s _1 is observed around 20:00 UT when 
the set of observations ends, which we identify in Section 
las Blob F. 

3.3. SOHO/LASCO white light observations 

SOHO /LASCO performs white light observations of 
the solar corona between 2.5 and 30 solar radii. Obser- 
vations of the 2000 June 28 CME were performed by 
both the C2 and C3 cameras for the duration of the 
event. The first C2 observation at 19:31:55 UT of the 
eruption showed a large plume of plasma off of the west 
limb (see Figure |4]). The flux rope observed by EIT 
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agation. The apparent velocity of this feature was ~300 
km s _1 at the low altitude end and ~600 km s _1 at the 
high altitude end while in the C2 field of view. 

The LASCO CME Catalog estimates a mass of ~7.3 x 
10 15 g and a kinetic energy of ~5.3 x 10 31 ergs. This cor- 
responds to a kinetic energy per unit mass of ~73 x 10 14 
ergs g _1 . However, the mass of this event is uncertain 
because this was a partial halo CME; furthermore, the 
kinetic energy is uncertain because it assumes that the 
entire plasma is propagating at the same velocity, and 
only considers the plane-of-sky velocity. 



Fig. 4. — LASCO C2 white light observations of this event, 
rescaled to emphasize faint structure. The approximate locations 
of the features observed by UVCS and analyzed using a time- 
dependent ionization code are labeled. The plane of sky velocities 
for these features were estimated using either O VI pumping infor- 
mation or complementary MLSO/MK4 observations. The precur- 
sor active region is denoted by the red plus sign. 



and MK4 had by this time broken up into several dif- 
ferent outwardly propagating blobs. In subsequent C2 
observations at 19:54:41 and 20:06:05 UT, each of the 
clouds propagated approximately radially outward from 
the location of the flare site and rising prominence with 
some expansion. The C2 image at 20:30:05 UT shows the 
loop identified as Blob F in Section |4] slowly rising above 
the center of the UVCS slit. This feature was probably 
ejected late during the eruption. 

Observations by both the C2 and C3 cameras show a 
bright stream of plasma propagating approximately ra- 
dially outward from the region of the flare site. Unfortu- 
nately, this feature passed just south of the UVCS slit. 
The feature exhibited significant expansion during prop- 



3.4. SOKO /UVCS observations 

SOHO /UVCS (jKohl et all fl995l . l2006h is a long slit 
spectrograph designed to study the solar corona at 1.5- 
10 solar radii. The use of both internal and external 
occulters keeps stray light at low enough levels to ob- 
serve the faint corona. The UVCS O VI channel con- 
tains two spectrometer light paths optimized for the O VI 
AA1032, 1038 doublet (primary) and Lya (redundant). 
Throughout the course of this event, the UVCS entrance 
slit was positioned at p = 2.32i? with a position angle of 
295° over the precursor active region near the northwest 
limb. The exposure time was two minutes, thus allow- 
ing high time resolution observations of the ejecta. The 
slit width was 49 /am. The observations are available in 
three wavelength ranges (denoted panels): 976-979 A, 
which contains the C ill A977 line; 1024-1045 A, which 
contains the O VI A1032,1038 doublet and H I Ly/3; and 
1210-1220 A, which contains H I Lya and the forbid- 
den and inter combination lines of O V at 1213 A and 
1218 A, respectively. Because the positions on the detec- 
tor overlap between the different channels, the primary 
channel N ill A989,991 doublet appears in the Lya panel 
at different wavelengths depending on the redshift. Be- 
cause of telemetry limitations, the UVCS observations 
were binned in three pixels in the spatial direction for a 
resolution of 21" per bin, and two pixels in the spectral 
direction for a resolution of 0.198 A (57 km s _1 ) for the 
C ill and O VI panels and 0.183 A (45 km s" 1 ) in the 
Lya channel. 

The UVC S observations of this event are summarized 
in detail by iCiaravella et al.l ( I2C 



Starting in the ex- 
posure at 18:59, UVCS detected faint, diffuse, and broad 
blueshifted O VI emission at position angles of 273-295°, 
indicating passage of the CME front. Knots bright in 
O VI and Lya first appear at 19:06 UT at a position angle 
of 273°. At 19:12 UT, a blueshifted bright knot in O VI 
and Lya appears at 277° (see Figure [5j this is associated 
with blob A) . This indicates the passage of cooler promi- 
nence plasma. By 19:14 UT, these knots fill the spatial 
direction along the slit between P. A. = 272° and 290°. 
Shortly after, two knots persist in the UVCS observations 
but the southern knot fades by about 19:30 UT. For 
many of the features early in the event, Ly/3, [O v], and 
C ill were off-panel because of the very large blueshifts. 
The northern knot persists for a while but drifts south- 
ward. Between 20:09 and 20:13 UT, a new diagonal or 
shear flow feature appears between P. A. = 288°-297° 
that has apparent O v], O VI, and C ill emission but lit- 
tle Lya or Ly/3. A few plasma streams persist afterward 
as the corona settles back into an undisturbed state. 
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Fig. 5. — UVCS observations of each of the features identified 
in Section [4] rescaled to emphasize faint structure. The pre-CME 
average background has been subtracted. From left to right are 
the C ill panel, the O vi panel, and the Lya panel. The rest 
wavelengths of the features are shown in the top row. 



3.5. GOES X-ray observations 

A C class flare was observed by GOES starting at 
18:48 UT with the decay continuing until about two 
hours after the event. Using the GOES 1-8 A light 
curve and subtracting the background flux, we estimate 
that the amount of flare energy emitted in this band 
is Ex ~ 3 x 10 28 ergs over the entire eve nt. Assum- 
ing L tot /Lx = 100 (e.g., lEmslie et al.ll2005[ ), this corre- 
sponds to an upper limit of flare energy averaged over 
the CME of <4 x 10 14 ergs g" 1 . This quantity is about 
an order of magnitude less than the kinetic energy per 
unit mass derived from the LASCO CME catalog. 

3.6. Yohkoh/SXT observations 

Yohkoh /SXT (jTsuneta et all Il991f ) performed pre- 
flare observations including the precursor active region 
(AR 9046) from 12:00 UT until 13:54 UT and from 
15:26 UT until 16:35 UT. Post-flare observations were 
taken from 19:26-19:46 UT. All three data sets used the 
thin aluminum filter. There were no other bright X-ray 
sources in the SXT field of view during this time. No ob- 
servations were taken during the peak intensity phase of 
the GOES flare. The observation at 16:50 UT shows that 
much of the soft X-ray emission from the precursor active 
region was behind the limb. The observations at 19:26 
and 19:50 UT, however, show some new bright X-ray 
emission in front of the limb. Because of the gap in SXT 
observations, we cannot rule out the possibility that the 
flare site was behind the limb and th us partially occulted . 
However, as shown most recently bv lAarnio et al.l (|2010), 
it is possible to have powerful and mas sive CMEs associ- 
ated with relatively weak flares (see also lReeves fc Forbesl 
2005D. 



3.7. SORO/MDI observations 

SOHO/MDl (jScherrer et al.lll995f ) takes high spectral 
resolution images of the Ni I A6768 absorption line to 
characterize velocity oscillations and Zeeman splitting in 
the photosphere. Full disk magnetograms were taken ev- 
ery 96 minutes. A subsection of the magnetogram was 
taken to match the field of view and orientation of the 
SXT observations. We use level 1.8 line-of-sight mag- 
netograms to measure the magnetic field and then esti- 
mate the magnetic energy of the active region prior to the 
event. We have included in our analysis magnetograms 
for 2000 June 25-28. On 2000 June 28, the flaring active 
region was on the limb. Due to the lack of resolution near 
the limb and in order to better characterize the magnetic 
energy associated with the active region, we include in 
our analysis magnetograms from 2000 June 25 when AR 
9046 was on disk. 

4. FEATURE IDENTIFICATION AND 
CHARACTERISTICS 

In this section we identify six features observed by both 
UVCS and LASCO with good density diagnostics. The 
positions and UVCS observation times of these features 
are shown in Table [2j along with the lines which were 
detected, undetected, or off of the UVCS panels. Ad- 
ditional properties of these features are shown in Table 

m 

Blobs A and B correspond to the first clear detections 
of CME prominence plasma at their respective positions 
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TABLE 2 
Features observed by UVCS 



Blob Time P.A. p/R Q 



Observed lines 



Upper limits 



Off-panel 



A 19:12 276° 2.45 O vi, Lya, N ill 

B 19:14 287° 2.35 O VI, Lya, C II, N in A992 

C 19:16 283° 2.37 O vi, Lya, C II 

D 19:21 284° 2.36 O VI, Lya 

E 19:27 285° 2.35 O VI, O v], [O v], Lya, C II A1037.0 

F 20:11 294° 2.32 O vi, O v], [O v], C in, N in A992 



C II, O v] [O v], C in, Ly(3 

O v], N in A990 [O v], C III, Ly/3 

O v], N in [O v], C in, Ly/3 

O v], C II, N in [O v], C in, Ly/3 
C II A1036.3, N in C in, Ly/3 

Lya, Ly/3, C II, N III A990 



TABLE 3 
Blobs Observed by UVCS 









O VI 






O V 




Lya 


LASCO 


Blob No. 


^1032/^1037 r max 


Vlos 


Vpos 


log Tig 


^1213/^1218 


log n e 


^max 


N e,max 


A 


1.72 


30 


-759 


1465 c 


7.0 a 






2.0 


23 


B 


2.21 


24 


-730 


1656 b 


7.5 a 






2.5 


19 


C 


2.60 


39 


-735 


1654 b 


7.0 a 






2.5 


22 


D 


1.67 


29 


-616 


1531 c 


6.9 a 






4.3 


20 


E 


3.73 


8.7 


-470 


1748 b 


6.6 a 


0.16-0.5 


6.2-6.8 


0.9 


53 


F 


2.14 


20 


-331 


250 




0.16 


6.8 a 




7.5 


Note. 


Key quantities 


derived for each of the 


features 


presented in 


Section [4] 


^1032 


//1037 and 



-^1213/^1218 are the line ratios used for the O VI and O v density diagnostics, respectively, described in 
Section[2] The line-of-sight and plane-of-sky velocities, Vlos and Vpos, are given in km s _1 . The plane-of- 
sky velocities were derived using O VI pumping information for blobs A E and MK4 observations for blob 
F. The number densities derived from these diagnostics are in units of cm -3 . T max is the maximum tem- 
perature in MK for each species allowed from UVCS line widths. -/V e;max is the maximum allowed column 
density from LASCO observations in units of 10 16 cm -2 . 
a Adopted density. 

b Indicates pumping of O VI A1031.91 by Ly/3 at velocities near 1810 km s _1 . 

c Indicates pumping of O vi A1037.61 by O VI A1031.91 at velocities near 1650 km s -1 . 



along the UVCS slit. There is strong O VI and Lya 
emission, with probable detections of N ill. The large 
blueshifts of 759 and 730 km s _1 , respectively, indicate 
that Ly/3, C ill, and [O v] are off of the UVCS panel. 
The O v] line is obscured by bright Lya, but upper limits 
of O v] are obtained. The densities are found through 
radiative pumping of the O VI doublet to be ~1 x 10 7 
cm -3 for blob A and ~3 x 10 7 cm -3 for blob B. 

Blobs C and D are located spatially between the lo- 
cations of blobs A and B along the UVCS slit but are 
observed by UVCS a few minutes later. The Ly/3, [O v], 
and C III lines are once again off-panel but upper lim- 
its of O v] are still possible despite bright Lya. Again, 
radiative pumping of O VI is used to estimate number 
densities of ~1 x 10 7 cm -3 for blob C and ~8 x 10 6 
cm -3 for blob D. 

Blob E is the first feature observed by UVCS with clear 
detections of both [O v] and O v] , allowing the density 
to be found using the intensity ratio between these two 
lines. The ratio of the O VI lines is iio32/^i038 = 3.73, 
which is substantially different from the collisional ratio 
of 2:1 and is indicative of significant radiative pumping 
of the O VI A1032 line. Pumping of this line by chro- 
mospheric Ly/3 indicates a total velocity of about 1810 
km s _1 , which is large compared to Flos ~ —470 km 
s _1 . The density inferred from O VI radiative pumping 
is n e ~ 4x 10 6 cm -3 . This is within the range of densities 
derived from the O V line ratio of 1.5-6 x 10 6 cm -3 . This 
interpretation implies that Vpos ^> Vlos, that this ma- 
terial is propagating close to the plane of sky, and that 
this plasma was ejected perhaps ~ 10-20 minutes after 
the initial eruption. 



A less likely possibility to explain the O VI ratio ob- 
served in Blob E is that the O VI A1032 line was being 
radiatively pumped by chromospheric O I AA1027, 1028 
emission at velocities near 1100 and 1300 km s _1 . Ac- 
cording to the SUMER sp e ctral atlas of the solar disk 
presented by iCurdt et al.l (j200l[ ). these two lines are 
about an order of magnitude fainter than chromospheric 
Ly/3. Consequently, this interpretation would imply that 
the number density is an order of magnitude smaller than 
inferred using the assumption of radiative pumping of 
O VI A 1032 by Ly/3. Because the number density de- 
rived using the assumption of radiative pumping by Ly/3 
is consistent with the number density derived from the 
O V line ratio, we conclude that radiative pumping by 
Ly/3 is much more likely. 

Blob F is observed by UVCS much later in the event 
(20:11 UT). We identify Blob F as the rising filament 
structure observed earlier by EIT and MK4. This blob 
appears as a diagonal shear flow feature in UVCS in O V 
and O VI emission. There is C ill emission with a differ- 
ent morphology, suggesting that the cool gas is included 
at least partiall y in a different component along th e line 
of sight (see also lAkmal et al.lfeOOlULee et al.ll2009f ). Lya 
and Ly/3 emission are both largely absent. Consequently, 
both [O v] and O v] are easily identified. The O V line 
ratio is used to derive a density of n e ~ 6 x 10 6 cm -3 . 

Column densities for each of these features are found 
using LASCO C2 observations rather than MK4 data 
because of the better signal to noise near the UVCS slit 
position. The mass per pixel for each C2 observation 
is found using the C2J4ASSIMG function provided by A. 
Vourlidas in Solar Soft IDL. Because of the relatively low 
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cadence of C2 observations, the C2 position for each fea- 
ture is estimated by assuming that the plasma from each 
feature continues propagating on the direct path between 
the flare site and the position on the UVCS slit where it 
was observed. The plane-of-sky velocity is assumed to be 
constant and is found using blueshifts and O VI pump- 
ing information (for blobs A-E) or MK4 observations for 
(blob F). C2J4ASSIMG assumes that the ejecta are propa- 
gating in the plane of the sky so we use the ratio of Vpos 
and Vlos to correct for the dilution factor and Thomson 
scattering angle. We assume that the column density in 
C2 observations goes down as /~ 2 , where I is the appar- 
ent distance from the flare site. Because there will be 
some departure from a constant velocity, we choose the 
largest column density within an apparent radius of 0.1 
to O.3i?0, depending on the estimated distance travelled 
between the UVCS and C2 observations. The results 
are shown in Table [3] and will be used to constrain the 
ionization fraction of O VI in the following sections. 

5. METHOD 

To estimate plasma heating rates for the features de- 
scribed in Section HJ we use a one-dimensional time- 
dependent ionization code to track the ionization states 
of the plasma between the flare site and when the fea- 
tures were observed by UVCS. The numeri c al me thod 
has been describ ed in detail bv lAkmal et al.l ([2001) and 
iLee et all ([2009), and we discuss key features here. 

Because we do not know the initial state of the plasma 
before the eruption, we run a grid of models with differ- 
ent initial densities, initial temperatures, heating rates, 
and heating par ameterizat ions. The initial densities are 
assumed to be in the range log(n/cm- 3 ) G [8.6,11.6]. 
The initial temperatures are assumed to be in the range 
log(T/K) G [4.6,6.6]. Consequently, the ejecta are al- 
lowed to be cool prominence plasma at relatively high 
densities or hot plasma from the ambient corona at lower 
densities. The final densities are presented in Tableland 
were found using one or both of the diagnostics discussed 
in Section [2] The determination of final density provides 
a significant constraint on parameter space. The velocity 
curve is scaled from the promine nce velocity curve fo r 
this event shown in Figure 2b of Ma ricic et al.l ((2006). 
This velocity curve includes an acceleration phase at low 
velocity. The ejecta are assumed to be expanding homol- 
ogously. 

We use four different heating parameterizations which 
we chose for simplicity and consistency with previous 
work. These simple parameterizations allow us to explore 
parameter space and the spatial dependence of heating. 
The first parameterization is the wave heating mod el for 
the fast solar wind presented bv lAllen et al.l (jl998). de- 
noted by Qahh oc exp (— d/O.ti?©), where d is height 
above the limb. This is physically motivated by wave 
heating models of the solar corona. While there may 
be departures from an exponential if, for example, the 
Alfven speed changes with altitude, this form is a rea- 
sonable approximation that allows for a gradual decrease 
of heating with height that does not depend on density 
or lead to excessive heating far from the flare site. The 
second parameterization is heating proportional to num- 
ber density, Q oc n, which also gives strong heating at 
low heights with a gradual decrease at higher altitudes. 
The third parameterization is heating proportional to the 



number density squared, Q oc n 2 , which concentrates 
heating at low heights where the density is large and has 
the same density dependence as radiative cooling. No 
physical mechanism is assumed for Q oc n and Q oc n 2 . A 
compilation of he ating models for coro nal loops is avail- 
able in Table 5 of lMandrini et all (|2000f ) with some show- 
ing heating proportional to a power of density, but it is 
not clear how applicable these models are to CMEs. The 
fourth parameterization is heating proportional to the 
time derivative of the sum of the kinetic and gravita- 
tional energies, Qkr oc c\(Uke + Uq) /dt. This expres- 
sion arises from Eq. 71 of KR when we assume for sim- 
plicity that a constant fraction of the released magnetic 
energy goes into heating the plasma within each model 
(e.g., dQ = —hdU m where h is the fraction of lost mag- 
netic energy appearing as heat which we assume to be 
a model parameter that is constant within each run and 
U m is the magnetic energy) Q 

For each model the ionization states are evolved using 
the relation 

dn z . _ N 

— = n e n z - 1 q i (Z,z - 1,T) 

-n e n z [ Qi (Z,z,T) + a r (Z,z,T)} 
+n e n z + 1 a r (Z,z + l,T), (4) 

where qi(Z, z,T) and a r (Z, z,T) are the ionization and 
recombination rate coefficients for an ion z of element Z 
at a temperature T. Initially, each model is in ionization 
equilibrium for the assumed starting temperature with 
coronal abundances. The particle distribution functions 
are assumed to be Maxwellian. The elements considered 
in this analysis are H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, 
and Fe. Cooling by radiative losses and adiabatic expan- 
sion are included in the analysis. A temperature floor is 
maintained at 5000 K, and model runs are rejected when 
the temperature exceeds 10 7 K. 

Once the grid of models is completed for each feature, 
the predicted line intensities are compared to UVCS ob- 
servations to determine which sets of parameters are ac- 
ceptable. The emission is assumed to come from the 
features bright in O V and O VI, while all other lines are 
used as upper limits since they might be emitted from a 
different component of the plasma along the line of sight. 
The UVCS observations alone constrain the line ratios, 
but the LASCO column densities constrain the ionization 
fraction of O VI, thus allowing us to place limits on the 
absolute line strengths. Line widths for O VI and Lya 
provide (usually uninteresting) upper limits on the final 
temperature. The results of this analysis are discussed 
in the following section. 

6. CONSTRAINTS ON PLASMA HEATING 
DERIVED FROM UVCS OBSERVATIONS 

Constraints on plasma heating for the features de- 
scribed in Section |4] derived using the techniques de- 

1 As pointed out by the referee of this paper, Eq. 72a of KR 
with (3 G ^ 0.22 implies that Eq. 72b should be 0.14 < h < 0.78. 
This lower limit of 14% of the magnetic energy available for heat- 
ing is in contrast to the lower limit of 58% reported by KR. 
We were able to reproduce Eq. 72a of KR using the expressions 
dUxE /&Um = — 2(1 — s) and (3q = —dUc/dUm- The latter expres- 
sion suggests that (3q and therefore h will in general be functions 
of time except when dU m /dt oc — dUo/dt (or, if we assume a point 
mass expression for gravitational energy, dUm/dt oc —V/R 2 ). 
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scribed in Section [5] are presented in Table 2J This table 
shows the components of the energy budget in terms of 
energy per unit mass, including the kinetic energy, gravi- 
tational potential energy, cumulative heating energy, and 
thermal energy of the plasma when observed by UVCS. 

The different heating parameterizations typically each 
have a characteristic temperature history pattern, as 
shown in Figure [6l For example, the allowed wave heat- 
ing models (Qahh) typically show the temperature drop- 
ping to <10 5 K before gradually heating up to a few 
times 10 5 K at the time when observed by UVCS (Fig- 
ure^). Heating proportional to number density (Q oc n) 
generally leads to a steady temperature after the ejecta 
leave the vicinity of the flare site (Figure [6}d). Occasion- 
ally, the temperature does approach the floor value for 
this parameterization during the middle of a run also. 
When Q oc n 2 , the ejecta tend to be heated to up to a 
few times 10 6 K before gradually adiabatically cooling 
to temperatures below 10 5 K (Figure [6b). The heating 
parameterization by KR does not greatly constrain the 
temperature history at early times but eventually tends 
to result in a gradually decreasing temperature at later 
times where heating does not provide enough energy to 
completely counter adiabatic cooling and radiative losses 
(Figure [6^). 

There are interesting upper limits on cumulative 
plasma heating for all features, and interesting lower lim- 
its on heating for blobs A and F (and to a lesser extent, 
blob E). For blobs A and E, the cumulative plasma heat- 
ing is constrained to be less than the inferred kinetic 
energy of each of the features. For blobs B-D, the cumu- 
lative plasma heating is constrained to be less than ~2-3 
the inferred kinetic energy. For blob F, the slowest fea- 
ture, the plasma heating is constrained to be comparable 
to or greater than the inferred kinetic energy. 

The results for blob E, the fast feature observed slightly 
later during the event, indicate that significant plasma 
heating comparable to the kinetic energy could only have 
occured early in the event when the ejected plasma was 
not far from the flare site. Heating proportional to the 
square of the density does allow heating up to 109 x 10 14 
ergs g _1 , but for this mechanism the heating is strongly 
concentrated at low heights where the density is high. 
The other parameterizations which allow for more grad- 
ual heating do not allow the cumulative heating to be 
greater than ~0.2 of the kinetic energy for this feature. 

Blob F, the slow feature observed late in the event, 
has the best constraints because it is the only feature 
for which the observed C ill line was completely on the 
UVCS panel. For Q ex Qahh, the temperature drops be- 
low 5 x 10 4 K for all of the allowed runs (Figure [7K). This 
behavior is possib le but unlikely, and was ruled out by 
lLandi et al.l (j201Qh for a separate event. For Q oc n, the 
allowed runs show steady or slowly increasing tempera- 
ture histories at a few times 10 5 K (Figure Et>). The ini- 
tial densities are 4 x 10 8 cm -3 , which is a coronal rather 
than prominence density at the low end of our assumed 
density distribution. However, the initial temperature 
is constrained to be <10 5 K for this parameterization, 
which would be unusual for plasma at coronal densities. 
The allowed models with Q oc Qkr show a gradually 
decreasing temperature history with initial densities of 
10 10 cm -3 , a broad range of initial temperatures, and a 
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Fig. 6. — Characteristic temperature histories for each of the 
different heating mechanisms shown using selected models for blob 
A. 

final temperature of 2 x 10 5 K (Figure [7h). 

7. ACTIVE REGION MAGNETIC ENERGY 

We estimate the total magnetic energy of the pre- and 
post-CME active region (AR 9046) using SOHO/MDl 
observations. We use an empirical relationship between 
the total unsigned magnetic flux and the magnetic en- 
ergy, 



2 x 10 32 



10 22 Mx 



1.35 



ergs, 



(5) 



presented by iFisher et al.l (|1998l ). The c oefficient for 
Equat ion [5] is estimated using Figure 8b of IFisher et al.l 
(1998). This method assumes that the magnetic field 
is potential and the results should thus be considered 
approximate lower limits or order of magnitude calcu- 
lations. It should be mentioned that the magnetic free 
energy is not accounted for but is probably of the same 
order of magnitude as the magnetic field energy found 
by assuming a potential field. 

Using the data set of magnetograms for 2000 June 25 
and Equation [5] we obtain an average magnetic energy of 
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TABLE 4 

Energy Budgets for Multiple Heating Models in Units of 10 14 ergs g~ 









Q oc ( 


3ahh 


Q 


oc n 


Q oc n 2 


Q oc 




Blob No. 


K.E. 


G.E. 


H.E. 


T.E. 


H.E. 


T.E. 


H.E. 


T.E. 


H.E. 


T.E. 


A 


136 (>29) 


7.4-7.8 


5.5-34.5 


0.62-4.1 


7.2-45.6 


0.31-2.0 


22.3-42.0 


0.03-0.05 


7.4-127 


0.10-1.6 


B 


164 (>27) 


7.9-8.1 


0.26-36.6 


0.03-4.8 


1.4-85.8 


0.07-4.25 


18.4-117 


0.03-0.21 


6.6-379 


0.1-4.6 


C 


164 (>27) 


7.7-8.1 


0.15-35.5 


0.02-4.6 


0.55-87 


0.03-4.3 


11.9-112 


0.02-0.2 


1.3-392 


0.02-4.7 


D 


136 (>19) 


7.9-8.1 


0.17-60.5 


0.02-7.5 


0.41-163 


0.02-7.4 


13.3-112 


0.02-0.15 


1.4-422 


0.02-5.6 


E 


164 (>11) 


8.2-8.3 


1.6-12.6 


0.2-1.6 


3.3-13.1 


0.16-0.64 


17.4-109.4 


0.03-0.2 


5.9-29.8 


0.09-0.44 


F 


8.6 (>5.5) 


5.5-8.2 


6.5-8.2 


0.67-0.95 


16.9 


0.9 






56.6 


0.41 


Note. - 


Components of the energy budget for the features observed by UVCS during this event, 


including the kinetic energy (K.E.), 



gravitational potential energy (G.E.), the cumulative heating energy (H.E.), and the thermal energy, shown for each of the heating 
parameterizations described in Section [5] assuming 10% helium. The kinetic energy is estimated using the velocities of O VI pumping 
for blobs A— E and MK4 observations for blob F, with the lower limits on the kinetic energies for all features found using VLos- The 
gravitational potential energy is given by GM@/R where R ranges from the apparent radius in the plane of the sky and the deprojected 
radius. 
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Fig. 7. — Temperature histories for the allowed models for blob 
F. 

8.6 x 10 31 ergs. Averaged over the day on 2000 June 28, 
the magnetic energy is found to be 1.8 x 10 31 ergs but 
is very uncertain because of the active region's proxim- 
ity to the limb. The magnetic energy estimates immedi- 
ately before and after the event do not significantly vary; 
this is expected because magnetic free energy is used to 
power the CME, and the CME itself should not signifi- 
cantly modify the line-tied magnetic field near the pho- 
tosphere. Using the LASCO CME Catalog's (uncertain) 
representative mass with the 2000 June 25 magnetic en- 
ergy estimate, this corresponds to ~110 x 10 14 ergs g _1 
and is comparable to the representative kinetic energy of 
the event. 

8. CANDIDATE HEATING MECHANISMS 



The most probable source of plasma heating during 
CMEs is the magnetic field. However, the mechanism 
which dissipates magnetic energy into thermal energy is 
not understood. In this section, we discuss several candi- 
date mechanisms which could be responsible for heating 
the ejected plasma. While quantitative model predic- 
tions for these mechanisms are generally not yet avail- 
able, we can provide observational constraints on some 
of these proposed mechanisms for this event. 

8.1. Heating by the CME current sheet 

Flux rope models of CMEs predict the formation of 
an elongated current sheet in the wake of the risin g plas - 
moid (e.g. JKopp fc Pneumaril 197ft iLin fe Forbesll2000[ ). 
Features identified as current sheets h ave been obs erved 
during, a number of CMEs (e . g . . iCiaravella et all 
Ciaravella & Ravmond 2008; Schettino et al.l 



2002 



2010; 



2011 



Savage et al.l 120101: iPatsourakos fc Vourlidasl 
Reeves fc Golubl l2011f ). These current sheets are ex 



pected to be unstable to the formatio n of plasmoids 
which may facilitate fast reconnection ([L oureiro et al.l 
20071: iSamtanev et al.l 120091: iBhattachariee et al.l 120091: 
Huan g & Bhattachari ej 12010b IShepherd fc Ca ssak 
2010; Ni et al. 201(1 lUzdenskv et alTl2Ql(3 : iBarta et al.l 
2008, l2QlQh . These current sheets have the potential 
to increase the thermal energy content of CME s via 
antisunward-directed exha ust. Recent anal ytical dSeatonl 
2008; IMurphv et al. [2010) and numerical (Re eves et al.l 
2010; IMurphv! l2010f ) studies predict that most of the 
mass, momentum, and energy flux from the CME 
current she et will be d i rected upward towards the rising 
plasmoid. ILin et all (j2004f ) argue that current sheet 
exhaust piles up in and around the rising flux rope so 
that the final parameters for CME evolution are set at 
heights of several solar radii. 

Outflow from the CME current sheet has the poten- 
tial to account for two key observed features of CMEs: 
that CMEs continue to be heated far from the flare 
site, and that the masses of CMEs tend to increase 
with time at distances of up to ~3-l OR ^ fro m the Sun 
(jBemporad et al.ll2007t IVourlidas et al.ll2010D . However, 
current observations do not strongly constrain the mass 
or energy budgets of CME current sheets. In particular, 
it is not known whether CME current sheets contribute 
substantially to the mass and energy budgets of the CME 
as a whole as some simulations and theories suggest (e.g., 
ILin et al.ll200i iReeves et al.ll2010D . 
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The predictions of the CME current sheet heating 
mechanism depend strongly on the i mportance of mix - 
ing and transport. In the model by [Lin et all ((2004), 
transport of reconnection outflow around the flux rope 
is assumed to fill each newly reconnected flux surface 
quickly without mixing plasma from different flux sur- 
faces. However, the reconnection outflow jets could pen- 
etrate deep into the flux rope and lead to significant 
turbulence and mixing (J. Karpen, private communica- 
tion, 2010). Mixing acts to spread the effects of heating 
over a larger volume. If mixi ng effects are not impor- 
tant, then the lLin et al.l ([20041 ) model predicts reconnec- 
tion heated plasma surrounding a core of cool material. 
iShiota et al.l (|2005l ) present simulations which suggest 
that slow shocks can permeate the flux rope, thus further 
contributing to heating. Further quantitative predictions 
of this heating mechanism can be made using a combi- 
nation of numerical simulations and analyti c theory. 

Co ol cores are commonly observed (e.g., lAkmal et al.l 
l2001f h suggesting that the cool plasma may be magnet- 
ically isolated from hotter nearby plasma. Some events 
do show significant he ating in the core of the CME (e.g., 
iBemporad et al.l l2007). In situ meas urements of low ion- 
izatio n state plasma in ICMEs (e.g. jLepri fc Zurbuchenl 
l2010f ) also provide constraints on mixing below heights 
where ionization freezes in. However, observations of the 
2000 June 28 CME do not provide significant constraints 
on the efficacy of heating by the CME current sheet. 

Extending the analysis performed in this paper to a 
large number of events will provide some information 
on the spatial dependence of CME heating. However, 
none the heating parameterizations described in Section 
[5] directly correspond to heating of ejecta by the CME 
current sheet . Num erical simulations such as those by 
I Reeves et al.l ([20101 ) can be used to investigate the effi- 
cacy of this mechanism, especially when used in conjunc- 
tion with a time-dependent ionization code to facilitate 
a comparison to observations. 

8.2. The kink instability 

The kink instability is driven by current parallel to the 
magnetic field and causes long flux ropes to develop a 
characteristic helical shape. Kinking is frequently ob- 
served in prominences and during solar eruptions (e.g., 
iRust fc LaBontdl2005l ). L ine- tying has a stabilizin g ef- 
fect on the kink instability ([Huang et al.ll2006ll20Tot ). but 
this mode may be destabilized by flux rope curvature. 
This instability has the potential to heat CME plasma 
by injecting turbulence through large scale motions. The 
turbulence then dissipates as the CME propagates away 
from the Sun. If this mechanism is most important, then 
it is likely that heating would occur over a turbulent dis- 
sipation time scale and would be greater in the vicinity 
of the flux rope. 

LAS CO and MK4 observations show that the flux rope 
becomes twisted. This behavior is consistent with the 
kink instability, but could be due to different effects. 
However, the transverse motions associated with this 
kink instability are much slower than the outward prop- 
agation of the flux rope. Thus the kink instability, by 
itself, is not likely to release enough energy through mass 
motions to explain total heating comparable to or greater 
than the total kinetic energy of a blob, as observed for 
several features during this and other events. However, 



additional magnetic energy could be released through re- 
connection events driven by the kinking behavior. This 
process is analogous to the kink-tearing behavior dur- 
ing sawteeth in tokamaks and other laboratory plasma 
confinement devices. 

8.3. Small-scale magnetic reconnection 

Th e ca ndid ate h eating mechanisms described in Sec- 
tions 18.11 and 18.21 are intrinsically linked to large-scale 
CME dynamics. An alternative to these models is heat- 
ing through small-scale, three-dimensional reconnection 
events which might or might not be driven by global dy- 
namics. For exam ple, flux rope expan sion models such as 
those by KR and W ang et all (|2009l ) predict that a sig- 
nificant fraction of the magnetic energy is converted into 
kinetic and thermal energy. The mechanisms by which 
this process occurs are not specified for these models but 
are very likely to be some form of magnetic dissipation 
(likely throu gh some com bination of reconnection and 
turbulence). lOwensI (|2009l ) presents a model which shows 
how internal reconnection can o ccur as a natur al result 
of flux rope expansion fsee also lXu et al.l 1201 If ). Small- 
scale reconnection is probably needed for CMEs to relax 
from a complex structure to the Lundquist configurations 
often observed in interplantary magnetic clouds (e.g., 
iLvnch et al.l I2004J ). Heating by small-scale reconnec- 
tion is analogous to the nanoflare model of coronal heat- 
ing. One possible manifestation of th is candidate heat - 
ing mechanism is the tearing mode ([Furth et al.| [l963). 
The tearing mode occurs frequently in magnetically con- 
fined laboratory plasmas in a variety of configurations 
and can moreover be driven by the kink instability. As in 
the case of the kink instability, line-tying where the flux 
rope is attached to the photosphere provides a stabiliz- 
ing effect and can change the eigenm ode structure and 
how t his mode scales wi th resistivity ([Huang fc Zweibel 
2009; lHuang et al1l2010l ). Tearing behavior can give rise 
to turbulence which can behave as an effective hyperre- 
sistivity (e.g. . Ivan Ballegooijen fc Cran mer 2008; St raussl 
fl988h . 

Because these processes occur on small-scales and mag- 
netic fields during CMEs are very difficult to diagnose, 
there are few observational constraints for this candi- 
date mechanism. During the 2000 June 28 CME, the 
O VI widths in blob F are ^100 km s -1 . This corre- 
sponds to an upper limit on the turbulent energy density 
of ~0.5 x 10 14 ergs g _1 , although there is some contri- 
bution to the line width from thermal broadening and 
perhaps shear flow. This upper limit indicates that tur- 
bulence must be continually ejected into the system and 
dissipated on time scales much shorter than the CME 
propagation/expansion time scale. Observational signa- 
tures of this mechanism include Alfvenic outflows and 
heating concentrated near the regions of maximum shear 
in the flux rope. 

8.4. Wave heating 

One of the leading mechanisms for heating of the am- 
bient corona and active regions is the damping of MHD 
waves driven by photospheric motions. The exponential 
heating parameterization used in this analysis was de- 
rived in the context of wave heating. We find that heat- 
ing rates of >100 and >63 times the coronal hole heating 
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rate of lAllen et al.l (| 19981 ) are required for blobs A and F, 
respectively, to explain the observed emission with this 
heating mechanism. This is consistent w ith the model 
behavior observed by lLandi et~a71 (j2010l), who s howed 
that heating rates >1500 times the lAllen et all (|1998f ) 
heating rate were required with this model during the 
2008 April 9 CME. For blob F, each of the allowed wave 
heating models shows a drop in temperature <5 x 10 4 K, 
which is possible but unlikely. These results provide fur- 
ther evidence that dissipation of MHD waves driven by 
photospheric motions is not able to explain the heating 
of CME plasma. 

Alternatively, recent laboratory experiments of the 
eruption of arched magnetic flux ropes show that in- 
tense fast magnetosonic waves resulting from the erup- 
tion are capable of tranferring energy to and heat- 
i ng the plasma in and aro und the erupting flux rope 
([Tripathi fc Gekelmanll20To[ ) . This mechanism might oc- 
cur in CMEs, possibly in conjunction with a large scale 
instability of the flux rope. The observations of the 2000 
June 28 CME do not provide meaningful constraints for 
this mechanism, except for the upper limit on turbulent 
energy density at UVCS heights discussed in Section [831 

8.5. Thermal conduction 

Thermal conduction along magnetic field lines is quick 
and therefore a potential contributor to the heating of 
CME plasma. The efficacy of this mechanism is limited, 
however, because of the shor t CME propagatio n time 
scales and long length scales. Landi et all (j2Q10h exam- 
ine thermal conduction in the 2008 April 9 CME and 
find that heating due to conduction is not sufficient to 
explain the observed heating rates. 

As discussed in Section [33J the 2000 June 28 CME is 
associated with a relatively weak C class flare. Assum- 
ing that the CME is heated approximately uniformly and 
that the flare was not strongly occulted by the solar limb, 
the lower limits on plasma heating for blobs A and F sug- 
gest that essentially all of the flare energy available must 
go into CME heating for this mechanism to be impor- 
tant. Therefore, it is unlikely that thermal conduction 
between the flare site and the ejecta is responsible for the 
inferred plasma heating. 

8.6. Energetic particles 

Some contribution to CME heating may be due to en- 
ergetic particles that were accelerated during the impul- 
sive phase of the event. A clear understanding of this 
mechanism requires a detailed description of the deposi- 
tion of energy by energetic particles into the bulk plasma 
(see, for example. lAllred et al.ll2QQ5h . The efficacy of this 
mechanism may be limited because the particle acceler- 
ation phase of flares is short compared to CME propa- 
gation time scales, and we know that heating continues 
at large distances from the flare site for many events. 

If the flare associated with the 2000 June 28 CME was 
not behind the limb, then it would be reasonable to infer 
that relatively few energetic particles were accelerated 
during the impulsive phase of this event. Additionally, 
blob F was observed ~80 minutes after flare onset which 
is probably too late for energetic particles accelerated 
from the flare to do much heating. Thus we conclude 
that energetic particles from the flare site are not likely 
to be responsible for substantially heating the ejecta. 



One complicating factor of energetic particles for this 
form of analysis is that non-Maxwellian particle distribu- 
tions can increase the ionization rates substantially in ad- 
dition to heating the thermal component of the plasma. 
It would be interesting in future work to see how a non- 
Maxwellian tail in the particle distribution function could 
affect the inferred heating rates. 

8.7. Heating by counteracting flows 

Filip pov fc Koutchmvl (j2002f ) describe a heating mech- 
anism for erupting prominences where upward concave 
flux rope segments yield shocks from colliding flows accel- 
erated by gravity. Because of the dependence on gravity, 
this mechanism is limited by the amount of gravitational 
energy available unless the plasma drops and rises multi- 
ple times. Between 1.1 and 3.Oi? , there is a difference in 
gravitational potential energy of llxlO^ergsg- 1 . This 
is much less than the kinetic energies of most of the fea- 
tures but greater than the lower limits on plasma heating 
for all of the features. For blob F, this is true only for 
the wave heating mechanism. However, the efficiency of 
this process in terms of available gravitational energy is 
probably less than unity since the distance plasma falls 
is likely to be less than the total distance available, and 
only a fraction of the plasma in a flux rope is likely to 
fall. Thus, unless the ejected plasma rises and falls mul- 
tiple times, we conclude that this mechanism is unlikely 
to account for sufficient heating during this event. 

8.8. Ohmic heating from net current in the flux rope 

The model bv iChenl (|1996l ) suggests that the injection 
of magnetic flux into a flux rope is a trigger for CMEs and 
that there is a net current through the flux r ope. Assum - 
ing a characteristic current of ~10 n A (e.g.. IChenl fl996). 
a flux rope minor radius of ~0.25i?©, a propagation time 
of ~10 3 s, and an electrical diffusivity of n ~ 10 4 cm 2 s _1 , 
we derive an expected heating rate of ~10 5 ergs g _1 . It 
is not surprising that this estimate is nine orders of mag- 
nitude below the lower bounds on blobs A and F, and 
thus cannot explain the inferred heating. 

The above estimate makes two implicit assumptions. 
First, the current density is assumed to be roughly uni- 
form across the flux rope cross section. Alternatively, 
the current density could be very strongly concentrated, 
which would indicate small-scale reconnection phenom- 
ena (see Sect ion [8^3]) . Second, we assume classical Spitzer 
resistivity. An anomalous resistivity due to wave-particle 
interactions may be present, but would need to be many 
order s of magnitude larger than the Spitzer resistivity 
(e.g.. iRakowski et al.ll2011h . Anomalous resistivity gen- 
erally requires large currents or sharp gradients over 
short length scales (i.e., the ion inertial length or ion 
sound gyroradius). Again, strong current filamenta- 
tion would be necessary and this would indicate small- 
scale reconnection phenomena as a more relevant model. 
Moreover, it is not clear whether the processes which lead 
to anomalous resistivity would produce volumetric heat- 
ing best described as being proportional to the square of 
the the current density. For these reasons we conclude 
that resistive heating from net current in the flux rope is 
not a viable means of adequately heating CME plasma. 

9. CONCLUSIONS 



CME ENERGY BUDGET 
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In this paper we use a one-dimensional time-dependent 
ionization code to investigate the energy budget of a 
CME observed by SOHO /UVCS. By running a grid of 
models with different initial densities, initial tempera- 
tures, and heating parameterizations, we are able to con- 
strain the total amount of heat deposited into the ejected 
plasma to counter radiative losses and adiabatic cooling. 

We perform this analysis for six features observed by 
UVCS and LASCO. Number densities of the plasma 
observed by UVCS are found either through radia- 
tive pumping of the O VI AA1Q32, 1038 dou blet (e.g., 
iRavmond fc Ciaravellall2004t iNoci et al.lll987l ) or by the 
classical density s ensitive O V AA1213, 1218 doublet (e.g., 
lAkmal et al.ll200"lt ). Both of these diagnostics are avail- 
able for one feature (blob E) and the derived number 
densities are consistent to within the expected system- 
atic errors. Total velocity information is found using ra- 
diative pumping of the O VI doublet and through white 
light observations. 

For two of the features (blobs A and E), the cumulative 
plasma heating is constrained to be less than or compara- 
ble to the inferred kinetic energy of the feature when ob- 
served by UVCS. For three of the features (blobs B-D), 
the cumulative heating is constrained to be less than ~2- 
3 times the inferred kinetic energy. For the slow feature 
observed late in the event (blob F), the plasma heating is 
constrained to be comparable to or greater than the ki- 
netic energy. Lower limits from two of the features (blobs 
A and F) yield cumulative heating energies of >5 x 10 14 
ergs g _1 . 

Next we discuss and consider constraints on a variety of 
heating mechanisms. Upflow from the current sheet that 
forms in the wake behind the rising flux rope could con- 
tribute substantially to both the mass and energy bud- 
gets of CMEs, but the energetics of CME current sheets 
are not well constrained for this or other events. The 
kink instability is able to drive turbulence by twisting of 
the flux rope, but the observed motions do not contain 
enough energy to heat the plasma for this event. Sec- 
ondary reconnection or tearing behavior (perhaps driven 
by the kink instability) can drive turbulence which dissi- 
pates and heats the plasma. Wave he ating; can occur ei- 
ther through photospheric motions (cf. lLandi et al.l feOlO) 
or by wave s generated by the eruption of the flux rope 
its elf (e.g.. [TV iDathi fc Gekelmari 120101 ). As also shown 
bv lLandi et al.l (|2010). wave heating by photospheric mo- 
tions is unlikely to be important for CMEs since the re- 
quired heating rates are orders of magnitude larger than 
those inferred from the observed nonthermal mass mo- 
tions. Thermal conduction can bring in thermal energy 
from the flare site or perhaps from the ambient corona. 
Energetic particles could deposit energy into the ejecta, 
but also increase ionization rates which would complicate 
this analysis. Because this event was associated with a 
weak (C class) flare, we consider heating by thermal con- 
duction or energetic particles unlikely to be important; 
however, the flare might have been partially occulted 
by the solar limb. iFilippov fc Koutchmvl (j2QQ2h suggest 
that heating could be due to colliding flows which were 
accelerated by gravity, but we conclude this is unlikely for 
this event since lower limits on heating for several fea- 
tures are comparable or greater than the gravitational 
energy available (see also lLandi et al.l l2010). 



We estimate the magnetic energy of the precursor ac- 
tive region using SO HO /MDI observa tions and an em- 
pirical relationship by Fis her et al.l (jl998[ ). This estimate 
is uncertain because the active region is near the limb, 
but a few days before the event the magnetic energy is es- 
timated to be ~8.6 x 10 31 ergs. This is the same order of 
magnitude as the representative kinetic energy from the 
LASCO CME catalog. This estimate assumes a potential 
magnetic field and thus is probably an underestimate. 

There are likely to be selection effects associated with 
the sample of events studied with the technique used 
in this paper (lAkmal et al.ll20QlHCiaravella et alJl200H : 
iLee et all 120091 : lLandi et al.ll2010D . The UVCS density 
diagnostics most often used are spectral lines of O V 
and O VI which are most prevalent in plasmas at tem- 
peratures of order ~10 5 K and most useful at number 
densities between 10 6 and 10 T cm -3 . Additional diag- 
nostics available for CME cores in the lower corona in- 
clude density sensitive line ratios of O IV, Mg VII, and 
Fe Vlll and are accessible with ins truments such as the 
EUV Imaging Spectrometer (EIS; ICulhane et all I2007D 
on Hinode. Thus with these analyses we miss plasma at 
cooler or hotter temperatures. For this event, we do not 
see morphological features in MK4 or LASCO observa- 
tions which do not have UVCS counterparts, but there 
may be additional components along the line of sight. 

Thus far this form of analysis has been used to con- 
strain CME heating rates one event at a time. How- 
ever, the model assumptions made by this and previ- 
ous works differ slightly, thus complicating attempts to 
make a systematic or statistical analysis of the problem. 
In future work, we will perform a standardized time- 
dependent ionization analysis for features in a large num- 
ber of events observed by SOHO /WCS. This will allow 
us to make direct comparisons of plasma heating rates 
during different events and provide tighter constraints on 
several mechanisms for CME heating. High cadence ob- 
servations by the Atmospheric Imaging Assembly (AIA) 
on the Solar Dynamics Observatory (SDO) are also well- 
suited for this analysis when appropriate density diagnos- 
tics are available. However, more detailed model predic- 
tions will be needed before most of the candidate mech- 
anisms could be definitively ruled out. 

In particular there are several open questions pertain- 
ing to the energetics of CMEs. These include: (1) What 
physical mechanisms are most responsible for heating 
CME plasma? (2) Are CME current sheets energeti- 
cally important to CMEs as a whole? (3) How uniform 
is heating within a CME? (4) How does CME heating 
depend on global CME properties such as speed, mass, 
and magnetic field strength and configuration? (5) What 
are the roles of energetic particles? (6) How do the ki- 
netic energy and cumulative heating energy compare to 
a CME's magnetic en ergy? (7) How do CME and ICME 
flux ropes "relax" (cf. lTavlorlfl986h ? (8) Is magnetic he- 
licity conserved during the evolution of these systems? 
We will address these questions in future work. 
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